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1 Abstract 

^ ! Friedmans cosmological equations for the scale factor are analyzed for the Universe contain- 
ing dark energy. The parameter of the equation of state of the dark energy is treated as an 
arbitrary constant whose value lies within the interval w G [—1.5, —0.5], the limits of which 
are set by current observations. A unified analytic solution is obtained for the scale factor 
Q . as a function of physical and conformal time. We obtain approximated solutions for scale 
^ ! factor to an accuracy of better then 1%. This accuracy is better then measurement errors 
I of global density parameters and therefore is suitable for the approximated models of our 
O^' Universe. An analitic solution is obtained for the scale factor in ACDM cosmological model 
P ■ both in physical and conformal time, for the description of the evolution of the Universe 
from the epoch of matter domination up to the infinite future. 
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2 Introduction 

Since the discovery of the accelerated expansion of our Universe [1], several theories have been 
CS ■ put forward to explain the phenomenon: presence of a cosmological constant, modification 
^ ! of gravity on the largest scales of space-time, presence of new light fields (see, for example, 
^ I reviews [2, 3, 4, 5, 6, 7]) and a Universe dominated by Chaplygin gas [8]. The substance 
that produces an effective anti-gravity force has been called dark energy. 

Where dark energy is modelled as the manifestation of a new physical field, it can be 
characterized by the equation of state p = wp, where the parameter w can differ from — 1 
^ , and, in the general case, is a function of time. In a simple model, where the dark energy is 
^ [ a scalar field with a positive energy (quintessence), the parameter w follows the limitation 
w > —1, while the cosmological constant corresponds to w = —1. However the equation 
of state can also be strongly negative, w < —1; dark energy with such an equation of state 
is called phantom energy. Current cosmological observations do not exclude the possibility 
of a time- dependent w where at relatively large values of red shift z the equation of state 
corresponded to a quintessence with w > —1, while in a later epoch it was manifested as a 
phantom energy with w < -1 [Sl[ini[IIl[nil2I]. 

The basic function characterizing the global properties of our Universe is the scale factor. 
An exact equation for this scale factor as a function of physical and conformal time is 
desirable for the analysis of several problems. This paper is dedicated to the derivation of 
such a solution and its approximated forms. 

The constitution of our Universe is characterized by parameters of density of the in- 
dividual components: Qm and Qq are densities of non-relativistic matter and dark energy, 
respectively fT5]. In the particular case of w = — 1 we define Qq = Q\, where A means 
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cosmological constant. There are several other components in the Universe, however their 
contribution to the total density is less than 1% (at all stages following recombination) and 
so will be ignored throughout this paper. In other words, a unified analytic solution for the 
scale factor is obtained for all epochs of the Universe since the moment of matter domination. 

We start out by describing the exact solution for the scale factor in the Standard Cosmo- 
logical Model - ACDM. Such an equation for the scale factor as a function of physical time 
was arrived to in the textbook [22] . 

We derive the solution for the scale factor as a function of both physical and conformal 
time. In this we calculate the allowed intervals of conformal time, as well as evaluate points 
of conformal time crucial to the description of the evolution of the Universe. 

In the second section we study the approximated functions describing the scale factor 
in the early Universe and the later stage, corresponding to A -term domination. We also 
analyze the level of certainty in these limits. 

In the following sections we find an exact implicit equation relating conformal time and 
the scale factor in a model with an arbitrary constant parameter —0.5 > w > —1.5, in the 
equation of state of the dark energy, as well as critical moments of time for such a solution. 
Further, we discuss the approximated solutions. 



3 Exact equations for the scale factor in the Standard 
Cosmological Model 

Analyzing the background space-time, we shall restrict ourselves to the standard model of 
the expanding Universe, filled with non-relativistic matter and using the A-term as the source 
of acceleration of cosmological expansion. The metric of the standard model is that of the a 
space-fiat Universe |l5j: 

ds"^ = dt^ - a^{t)d^^. 

The scale factor a{t) is determined from the Friedman equation, which can be written in the 
form: 



a{t) 



Hi 
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(1) 



where Hq is the current value of the Hubble constant, a(to) = ao = 1 in this paper is 
understood as the current value of the scale factor and a dotted variable indicates a derivative 
with respect to physical time t. Furthermore, we work with a system of units corresponding 
to c = 1. In this work we use VLrn = 0.27 and VL\ = 0.73 as the current values of the density 
parameters [H]. Here we have chosen a set of parameters based on 7-year observations of 
the WMAP probe and recommended for use. It must be noted that there exist several other 
models corresponding to different values of the global cosmological parameters. However, 
for our analysis this difference is not a principal one and a recalculation of the fundamental 
results is easily accomplished for any choice of these global parameters. 

This standard model of our Universe has been revaluated time and again wince the 
discovery of its accelerated expansion. 

Our aim is to find an explicit form of the scale factor as a function of physical time t and 
conformal time f], defined as: 

dt 
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We produce these equations and determine the exphcit dependence of the scale factor on 
physical and conformal time, as well as derive several often used equations determining the 
age of our Universe, moment of transition from decelerated to accelerated expansion, etc. 

Friedmans first equation may be written in the form: 



f da{t)\' 



Hi 



a3(t) 



(2) 



a2(t) \ dt J 

This equation can be integrated to obtain the following relationship between time and 
the scale factor: 



Hnt 



dx 



(3) 



Here the upper bound of integration depends on the scale factor as 



Xa 



-a 



3/2 



The current age of the Universe can be determined from this integral by substituting 
a = 1 or X, = 



n— . The age of the Universe as a function of these global parameters 
TT), 
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The solution for the scale factor can be calculated [22\ from the explicit form of integral 



(ED: 



a{t) 



1/3 



sinh 



1 2/3 



(5) 



We now analyze the value of the scale factor a{t) at the two limiting cases. The first 
is the limit of small Hot << 1. In this case, the contribution of the A-term to the total 
density is infinitesimally small and the Universal expands at a decelerating rate, by the law 
corresponding to matter domination. Then, taking the Taylor expansion of we have: 



a{t) 
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[Hot) 



2/3 



(6) 



which corresponds to the scale factor in the stage of matter domination. 

We now calculate the moment of physical time at the beginning of matter domination. 
This is done by equation of the parameters of density of matter and relativistic particles 
(those with an equation of state p = 3p). The value of physical time at the moment of equal 
densities of matter and radiation is found from the standard condition [15], according to 
which the scale factor a^, at this moment is equal to the ratio of densities of radiation at the 
current point of time to the current density of matter: 



dm. 
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Here is the density parameter of relativistic radiation. According to WMAP [16] the 
radshift of the epoch is Zm = 3196. Then the moment of physical time corresponding to 
the epoch of matter domination is: 

f /O 3/2 

I'm / i'A \ "m ^ij^ 



to \^mj /l + 

In I 1= 



In the other hmit, i.e. for large values of Hot » 1, the contribution of the A-term to 
the total density becomes dominant, while that of matter significantly diminishes. Now we 
obtain a different dependence of the scale factor on time: 
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ait) = i^-^j exp [V^AHot) , (8) 

which corresponds to the scale factor of a Universe with a de Sitter expansion law. 

We now examine the switch from decelerated to accelerated expansion. To this end we 
analyze Friedmans second law. In terms of physical time the second equation reads: 



1 fdMt)Y _ 1^2 



(9) 



a{t) \ dt^ J 2 ° 

The moment of time t = when decelerated expansion changes into accelerated expan- 
sion is determined by the equation: 

^ = 0. (10) 

At this moment of time the value of the scale factor is: 



1/3 



oaH^I . (11) 



which corresponds to the red shift: 



The moment of time when the stage of deceleration changes into the stage of acceleration 
is determined according to: 

i/„(.= /^^L=. (13) 

We calculate the values of physical time at various epochs by for global parameters 
Qm = 0.27, Q\ = 0.73 Hq = 71 km/s/Mpc. The moment of matter domination corresponds 
to value of red shift z = 3196 and constitutes 98357 years from the beginning of expansion. 
The moment of change from deceleration to acceleration is =7.1 billion years, and the red 
shift at this point has a value of z^ = 0.82. Finally, the age of the Universe constitutes 13.75 
billion years, with the current value of red shift equal to z = 0. The Universe infinitely in 
the future corresponds to a physical time t = oo. 



4 



The dependence of the scale factor on physical time is not sufficient for several formulae. 
To be used in these formulae, it must be expressed as a function of conformal time. This is 
most easily obtained by rewriting Friedmans first equation in conformal time: 



/ da{r]) \- 



Hi 



a'^{ri) \ drj J 
The solution to this equation in integral form is: 



(14) 



HoV = [ ^ (15) 

This integral can be expressed in explicit form as an elliptical integral F{ip, k) of the first 
kind. For this, the integral is rearranged as follows: 

u 



/Q n1/3 

The upper limit of integration here is m = ( jy^ ) o,- The integral ( fT6|) is solved [23] as: 



\ 1 + (1 + V3)m 2 / 



(17) 



Elliptic integrals can be inverted using Jacobi elliptic functions to obtain an explicit 
formulation of the scale factor as a function of conformal time. Thus we use the elliptic 
cosine [21] to describe the scale factor function: 

air,) = ( ^ ^-cn{y.k) 

\^aJ V3il + cn{y,k))-il-cniy,k)) 



where the argument of the cosine is y = S^^'^Q^^nU'^Hori and its modulus is k = — — ~ 
0.97. 

Now we need to determine the moments of conformal time corresponding to the following 
epochs: the begin on expansion rjs, the beginning of matter domination rjm, the epoch of 
change from decelerated to accelerated expansion r^A, the current epoch tjq and the infinite 
future T^oo- Since conformal time does not have a direct physical meaning, the way physical 
time does, we will be calculating not r] itself, but the two quantities y and HoT]. It must be 
noted that the quantity y can be evaluated exactly, while HoT] is expressed through y and 
the two observationally determined constants (density parameters) ^a, ^m, that have a 
precision of a few percent. The precision of NqT] is thus constrained by that of the density 
parameters, and we state the numerical values of this quantity up to only two decimal points. 
It follows that our model is inadequate for describing the early Universe, i.e. this model gives 
an inexact law of evolution of the scale factor in the interval from the beginning of expansion 
to the moment of matter domination. However, after this moment our solutions get more 
accurate as the point of time in question moves closer to the current moment. Thus, while 
we so analyze evolution starting at the beginning of expansion, the solution is considered 
qualitative until the moment of matter domination. 
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At the beginning of expansion the scale factor has a value of zero, which is thus the lower 
bound for allowed values of conformal time. In other words, at the beginning of expansion 
y = 0, Hor]s = 0. 

The value of conformal time at the moment of equal densities of matter and radiation is 
found from the standard condition, as in the case of physical time. This gives us the value 
for the variable ym = 0.05 and the value of conformal time at this moment: 

HoVm = 6.2 ■ 10-2. (19) 

This moment of time we shall henceforth call the beginning of the stage of matter dom- 
ination. 

We are now left to calculate conformal times corresponding to the switch from deceler- 
ating to accelerating expansion, the current epoch and the infinite future. The moment of 
time corresponding to (ITT]) is: 

1 + 2^/3 _ 3I/2 
cn{yA,k)= ^^^1/3^31/2 - 

The solution to this equation is y\ = 2.27, which for the accepted values of the global 
parameters is equivalent to: 

HoVa = 2.82. (20) 

The current epoch is obtained by setting the scale factor to one. This leads to the 
solution: 

cn{yo, k) 



For arbitrary values of the density parameters this solution can be rewritten as y^ = 2.78. 
Alternatively, for the selected values of these global parameters we have: 

Hovo = 3.45. (21) 

The moment of time in the infinite future corresponds to the scale factor going to infinity 
(in the standard cosmological model the Universe expands infinitely). Correspondingly the 
first zero of the denominator in equation ( !T8|) determines the moment of conformal time 
representing the infinite future. This solution is: 

1 - 3^/2 
cn{y^,k) = Y^^- 

Again, the solution to this equation for arbitrary values of the density parameters is yoo = 
3.69 and for the selected values gives: 

HoVoo = 4.57. (22) 
Thus the complete interval for which y can be defined is: 

ye [0,3.69], 

which leads to the interval of allowed values of Horj being: 

Hove [0,4.57]. 

The allowed interval of values of conformal time is thus established. 
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4 Approximate solutions, error analysis 



In discussing approximated solutions to the scale factor, it must be immediately noted that 
these approximations are justified only for z « 10^, i.e. at or after the stage of matter 
domination. At earlier stages the Universe is dominated by radiation, a scalar field of 
inflaton and possibly yet-unknown forms of matter. Therefore, all our calculations for the 
real Universe as justified in the interval H^rj G [0.0122,4.57]. Nevertheless, within this 
interval of conformal time we may use either the exact or the approximated equations for 
the scale factor, as is discussed in this section. 

Given an exact solution to Friedmans equations for a scalar factor (in terms of conformal 
time) we find an approximation to the scale factor in the limit of small values of a. Small 
values of a correspond to HqT] « 1. The elliptic cosine upon Taylor expansion has the 
form( [21! Pg-46, f.H): 

cn(y. A;) = 1 - 1^ + (1 + 4P)|| - (1 + AAk' + 16fc^)|J + ... (23) 

Leaving terms up to the quadratic in y we obtain the approximate expression for the scale 
factor: 

m = % {Hovf , (24) 

which corresponds to the expression for the scale factor in the stage of matter domination 
in terms of conformal time. 

Now we analyze the asymptote of a{ri) as it approaches r/oo- For this we define y = yoo + ip 
and Taylor expand the elliptic cosine about yoo- We first use the formula from addition 
theorems of the elliptic function: 



cn{yoo, k)cn{ilj, k) - sn{yoo, k)dn{yoo, k)sn{ilj, k)dn{ilj, k) 
1 — k'^sn'^{yoo, k)sn^{ip, k) 



cnyyoo + ip, k) — - 27" rr~727~n~\ • 



Here, sn{y, k) is the elliptic sine, while dniy^ k) is the third elliptic function defined by Jacobi 

Using the properties of the Jacobi elliptic functions we determine the values of the three 
functions at the point ?/oo, In the calculations below we shall not explicitly write out the 



2 + a/3 

modulus k, and assume its value to be = — — ^ everywhere. 



cn(?/o 



sn{y„ 



dn{y^ 



1 - 

i + Vs' 

2^3 
2 + VS' 



2 

Leaving the highest order term of ip in the Taylor expansion of both the numerator and 
denominator in (18), we obtain the following expression for the scale factor: 
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Plugging in the expression for in terms of 7] we finally obtain: 



"^^^ = TT^T^ ^' ^^^^ 



The time variable can be redefined by introducing the new variable fj = rj — r/oo- Then 
the quantity fj appears negative, and the scale factor as a function of rj is: 

a{fi) = - ^ 



which corresponds to the known solution for the scale factor (as a function of conformal 
time) in a de Sitter metric [T5] . 

We now analyze errors that occur when evaluating the scale factor through formulae ( 12^ 
and ( 126|) and matching together these approximated solutions, for example, at a ~ 1. We 
start with the approximation (!24|) . Evaluation of errors of the approximate solutions is done 
by plugging in ( |2^ at the moments of time evaluated in (l20l) and (I2T1) . The error evaluation 
at the remaining points of time is trivial. At the point ( IT9l) the approximate and exact values 
of the scale factor are equal with an accuracy of greater than 1% while at the point f l22|) the 
ratio of the approximated scale factor to the exact value equals zero. 

We first evaluate the accuracy of extrapolating the lower-limit approximation f l24|) to 
higher values of conformal time: 

OA ^ 21/3 ^ 0.94, 



OA 



1/3 

^ 0.81. 



do ^ (2.12)2 
ao 4 V ^A / 

It is now seen that the approximated curve described by formula (12^ lies under the curve 
described by the exact formula. The deviation between the two constitutes 6 % at the point 
HoTjA (See FigHJbut increases to 19% at HqTIq. 

Calculating the accuracy of interpolation of the upper-limit approximation fl2^ into lower 
conformal times, we have: 

— = 2^/Vl-07 ^ 1.18 

OA 

1 /n \ 1/3 



, ^ . 1-06 
ao 0.68 V^A. 

Here the approximated curve described by formula (l26l) . lies above the curve described 
by the exact formula. At HqIIq the deviation is 6% while it increases to reach 18% at HqT]/,^. 

Since the error upon extrapolation from small values of the scale factor that that upon 
interpolation from bigger values have different signs upon comparing with the exact value, 
the errors add up. Therefore the error upon lacing together the scale factor using these two 
approximated formulae constitutes almost 25%. 

An even greater error appears upon comparing the first and second derivatives of the 
scale factor. The first derivative of the logarithm of the scale factor from the side of small 
vales is two times bigger than the analogous derivative from the side of greater values of the 
scale factor. A comparison of the second derivatives gives a ratio of values: 



0.54, 
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here prime denotes a derivative with respect to conformal time. Thus the approximated 
equations for the scale factor should only be allowed in in calculations that require accuracy 
no greater than 25%. 

The approximate value of the scale factor can be analyzed close to a = 1. However, 
approximating the value of a close to unity by power functions is inefficient, since a large 
number of terms of the Taylor expansion would be required to obtain an accuracy comparable 
to 1%. 

Therefore, we shall use another approximate solution for elliptic functions of the form: 

27r ^ g"-V2 



n=l 



+ q 



Where K is the complete elliptic integral of the first kind along modulus fc, the auxiliary 

function q = exp{—^^^^^),and the primed value is the complete elliptic integral of the first 
kind from the remaining value of the modulus. 

The first three terms of this expansion already allow for errors less than 0.5% across the 
entire interval of rj. Therefore, if the problem lies in finding an approximate solution for the 
scale factor with accuracy of the order of 1%, it is sufficient to use the sum of these first 
three terms of expansion. The value of the complete elliptic integrals is K = 2.768, and 
K' = 1.598, while q = 0.163. 

One can represent the scale factor as: 

a{7]) = ai——, (27) 



here 



V3-1 

ai = 



1/3 



/ X T 27ryg /I vry g Siry q"^ bny 

aniri) = 1 -rr— COS — — H 5- COS — — H F COS — — 

^" kK \l + q 2K l + g3 2K 1 + q^ 2K 

a.(r,) = 2- v^+^(^^cos— + ^cos— + ^cos — 
let introduce variable ( = y/Q^Hot] and we get the approximate form of scale factor: 

1- 0.8159 cos 0.5361C- 0.1541 cos 1.6083C- 0.0252 cos 2.6805C , , 

air?) = 0.263 : : -, (28) 

^" 0.2680 + 0.8159 COS0.5361C- 0.1541 cos 1.6083C + 0.0252 COS2.6805C 

This equation approximates exact equation f|T8|) with accuracy better then 1% over the 
whole interval of allowed values of rj. 

5 General Case of Parameter w 

Analysis of observational data indicates that parameter w is close to minus unity [w = — 1) 
tl6j|, [HI, [IS]- An interval of acceptable value of w changes significantly depending of the 
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analysis method, observational errors, and bias. Therefore the value of w can be different 
from -1. Analysis of cosmological models with w 7^ — 1 is requeired. 

We now look at the dependence of the scale factor on time for an arbitrary value of 
w G [—1.5,-0.5]. This interval of values is slightly larger than that set by observational 
limits, however, we use it for convenience. 

The equation for conformal time r/ as a function of the scale factor a can be obtained 
from the inverse function. The solution 77(a) is: 

a 

HoV = [ , "^"^ (29) 

5.1 Analitical solution for the scale factor 

Finally, we write out the analytic function describing the conformal time as a function of 
the scale factor: 

iv^ff„„=VSF(^.l;^;-§^a»H) (30) 

2 \b\w\ 2 b\w\ \Lm J 

Here F{a, P;^; z) is the Gauss hypergeometric function. This equation must be inverted 
to obtain the equation for the scale factor as a function of conformal time. However, the 
authors are unaware of special functions that could solve this problem. Therefore, where 
such a solution is required, it is determined by numerical methods. FigJH shows the form of 
the function 0(77) for values of the parameter w = —0.5, —1.0, —1.5. 

Now we find the asymptotic behaviour of the at small and large values of conformal time. 
It must be noted that the first of these asymptotes must correspond to a matter dominated 
Universe, and the latter to a Universe dominated by dark energy. 

The Gauss hypergeometric function can be Taylor expanded on the argument — ■ry^a^'"''; 

in the case where the absolute value of this argument is less than unity, this expansion 
converges. We note that: 

1 1 6\w\ + 1 



;0 = 1. 



6\w\ ' 2' 6\w\ ' 

This means that for small values of the scale factor and thus conformal time we once again 
obtain formula fl2^ : 

which describes the evolution of the scale factor in the stage of matter domination. 

In the second case, where a ^ 1, functional relations of hypergeometric functions may 
be used to obtain a series expansion of the function F {a, (3; 7; z) in terms of hypergeometric 
functions of the variable 1/z [2^ : 

T('y)T(B-a) 1 

F{a, /3; 7; z) = { {-zf F{a, a - 7 + 1; a - /3 + 1; -) + 

r(/3)r(7-a) z 

r(7)r(a-/3), ^ _^ _ 



p. .^ V ^)'^F(/3,/3-7 + l;/3-a + l; 

r(«)r(7-/3) z 



Using this expansion and the equivalence of hypergeometric functions to unity for an 
;ument approaching zero: 1/z = — 
the limiting case of a large scale factor: 



argument approaching zero: 1/z = —jj^a ^l'"' — )■ 0, we obtain the following expression in 
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2 1 l-3|m| 

HqT] = HqT]^ - —- -^=a 2 , (31) 

3\w\ - 1 



Here HqT]^ is determined by formula (IMI) . Meanwhile, the scale factor is now determined 
by the expression: 

a(^) = ^ (32) 

^^^^^^ ^A/^\i^o(^oo - V) 

It may also be noted than for tt; = — 1, the equation (!30l) reduces to: 

1 



a{ri) 



Although the set of exact results relating to the evolution of the scale factor got an arbi- 
trary value of parameter w appears limited in completeness, certain interesting formulations 
may be suggested for the implicit expression of the scale factor through conformal time. 

The hypergeometric function of the given values of parameters may be expressed in 
terms of an incomplete beta function. This function is used in mathematical statistics. 
Some minor standard reformulations also allow us to express the hypergeometric function in 
terms of Legendre functions. 

We produce the equation for the scale factor in a form convenient for the analysis of its 
evolution in the limit a >> 1: 

3|iu|-l 

2 1 

(33) 



p(\\ L. ^ L. ^"^ ] - 1 ^^Hq (r/oo - r/)^^ 

V2'2 elwl' 2 Q\wy VLqa^\^\ 



Let us consider the conformal time intervals. It is easily shown that this integral (129|) 
has a finite value for infinite limits under the condition 3|w| > 1. Otherwise, the value of the 
integral (12^ diverges at infinity. In other words, in the interval of values of \w\ of interest 
to us the integral (1291) is finite and thus so is the finite interval of values of rj. 

The value of conformal time for which the scale factor goes to infinity is: 

The complete interval of conformal time depends on the parameter w. For an upper limit 
of the parameter w = —0.5 the interval of change of conformal time is i^o^oo ^ [0, 7.75]. For 
the lower bound w = —1.5, this interval is HqTIoq G [0,4.20]. 

Let us estimate the conformal time at the moment i/o'?? of zero decelaration. 

This moment is: 



1 1 



6 1 I 



2^^™^°^^=(3|^|-1)V6H {n-J ^ (,6H ' 2' ^ • ^^^^ 
The graph of HQrjg as function of argument w is shown in FigJSl 
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5.2 Measuring the parameter 

The red shift at the moment of time when deceleration turned to acceleration, in this general 
case, corresponds to: 



1 . ( , 



\ l/3|u-| 

^. = --1= ( (3k|-l)^ -1, (36) 

Qjq y i Lyyi J 

Measurement of the parameter of the equation of state poses a tough but very important 
problem. We focus on one possibility of measuring \w\ from measurements of red shift z^. 
Measuring the value of red shift upon which decelerated expansion shifts to an accelerated 
one, it is possible to measure the parameter w. It is interesting to note that the dependence 
is not monotonous. We also note that the maximum of the curve in Fig. [3] corresponds 
to the solution to the equation: 



\6\W\ — 1/ \lm 



and consequently there exists a unique relationship between the parameter w and the max- 
imal value of Zq. 

To estimate the errors of w parameter we have to evaluate the derivative: 

- ' + ^ ln(3|^| - 1) - + In (37) 



d\w\ 3w^ V 3|u;| - 1 U 



It is possible to try to approximate the solution for a scale factor with arbitrary w to 
that of the standard model of the Universe. To this end we determine the function with an 
accuracy of: 

A = <^M-<^.^ = -^)^ (38) 
a(?7, w = —1) 

The graph of this function if shown in FigJSl It is observed that the approximation varies 
from the exact solution by over 10%. This means that the solution to the scale factor in the 
standard model of the Universe is not suitable to describe the evolution of a Universe filled 
with dark matter with an arbitrary parameter w. 



6 Conclusion 

In this paper we analyzed exact solutions for the scale factor to model a Universe filled with 
matter and dark energy. 

We first comprehensively analyzed exact solutions for the scale factor in the standard 
model of the Universe filled with dust-like matter and containing a A-term. Exact solutions 
are derived and presented in terms of physical and conformal time. The latter exact solution, 
in conformal time, is expressed in the form of Jacobi elliptic functions. We consider it 
to be particularly important, since the analysis of fluctuations of the gravitational field 
is conducted in terms of conformal time. We also studied approximated solutions in the 
standard model of the Universe. It was shown that approximated solutions in the form 
of power laws were not satisfactory. Their inaccuracy significantly exceeds the error in 
determining fundamental global parameters of our Universe. Thus, it is impossible to use 
these approximated solutions for a majority of problems. We found an expression for the 
exact solution as a sum of three trigonometric functions, which comprise a solution for the 
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scale factor with an accuracy better than 1% in the entire interval of allowed values for 
conformal time. Such an approximation may be used for all problems, since the accuracy of 
the approximation is higher than the rms of global cosmological parameters. 

We also found exact solutions for conformal time as a function of the scale factor for 
arbitrary values of the parameter w. A satisfactory approximation of this solution could not 
be found in terms of reasonably simple functions, although certain useful approximations 
are discussed in the paper. 

The paper also discusses the entire allowed interval of values of conformal time for a 
Universe filled with dark energy with varying values of the parameter w, as well as certain 
other critical points of time: the moment of the beginning of matter domination, the moment 
of change from decelerated to accelerated expansion and the moment which corresponds to 
the infinite value of the scale factor. 
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Figure 1: The diagram shows the curves describing the evolution of the scale factor. Values 
of the scale factor are placed along the vertical axis and HqI] grows along the horizontal axis. 
The exact value of the scale factor is represented by the solid line, the approximated formula 
fl26|) is graphed by the dotted line, and the dashed line shows the scale factor as evaluated 
from formula (^^. See text for details. 
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Figure 2: The diagram shows the dependence of i^o^oo on the parameter w. it is observed 
that when the parameter is in the interval —0.8 < tu < —0.5 the value of i^o^oo varies 
strongly with the parameter. After -0.8 up to the end of the allowed interval of values for 
the parameter, HotJoo depends on it weakly, with its value staying between 4.98 and 4.2. 
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Figure 3: This diagram shows the value of red shift Zq upon the change from decelerated to 
accelerated expansion. The argument of this function is the parameter of equation of state 
\w\. It is interesting to note the non-monotonous dependence of Zg on w. The maximum 
value of Zq corresponds approximately to w = — 0.9. . 
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zeta 

Figure 4: Here scale functors as a function of Horj are plotted for three values ofdark energy 
parameter w. Bold line corresponds w — —1, asterics correspond to - w — —0.5, and light 
line corresponds to - w = —1.5. 
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Figure 5: This figure represents the dependence of Horjq on the parameter w. 



19 



Figure 6: This figure represents the function A of the variable zeta = HqI]. The thick vertical 
line represents the moment of time corresponding to the beginning of accelerated expansion 
of the Universe, while the next vertical line represents the current moment of time. It is seen 
that at the point of change from deceleration to acceleration, the value of A already exceeds 
10%. This means that for several problems, a Universe dominated by dark matter cannot 
be approximated by the standard model of the Universe. Instead, suitable approximations 
must be considered for each case. 
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